Abstract 



Central limit theorem for fluctuations of linear eigenvalue 
statistics of large random graphs. Diluted regime. 

M. Shcherbina^ B. Tirozzi* 
'''Institute for Low Temperature Physics, Ukr. Ac. Sci 
O ' *Department of Physics, Rome University "La Sapienza" 

> 
O 

^. 

, We study the linear eigenvalue statistics of large random graphs in the regimes when 

' the mean number of edges for each vertex tends to infinity. We prove that for a rather 

wide class of test functions the fluctuations of linear eigenvalue statistics converges in 
' ^ ' ' distribution to a Gaussian random variable with zero mean and variance which coincides 

with "non gaussian" part of the Wigner ensemble variance. 

-3 . 

! 1 Introduction 

In this paper we study the spectral properties of ensembles of adjacency matrices of large ran- 
dom graphs. Following Erdos (see, e.g. |6]), we introduce the probability measure considering 
J> ' the set of all graphs with n vertices and set the weight of each graph G as 

(N 
in 



P{G) = {pn/n)<^\l-pn/np)-<''\ (1.1) 



where e{G) is the number of edges of G and Q < Pn ^ n. The set of n-vertices graphs with this 
measure (usually denoted by G(n,p„/n)) is one of the classes of the prime reference in the 
theory of random graphs. Most of the random graphs studies are devoted to the cases where 
Pn/n — >• 0, as n — )• c«. There are two major asymptotic regimes: p„ S> 1 and pn = 0(1) 
and corresponding models can be called dilute random graphs and sparse random graphs, 
^ ■ respectively. 

It is well known that there is one-to-one correspondence between the graphs and their 
adjacency matrices. For G{n,pn/n) the ensemble corresponding to (II. ip consists of random 
symmetric n x n adjacency matrices ^ is ^ = {aij}fj^^ with an = 0, and i.i.d. 

1, with probability Pn/n, ,^ 
0, with probability 1 — pn/n, 

This is a particular case of the random matrix ensemble. Since the pioneering works by "Wigner 
|19) a big part of the random matrix theory is devoted to the limiting transition n — )• oo. The 
results obtained with this limiting transition provide a rather good approximation of the 
spectral properties of random matrices (or random graphs) of a finite dimensionality. 

An important advantage of random matrices (jl.2p is that their entries are independent up 
to the symmetry condition (ojj = Uji). This allows one to use the methods of random matrix 
theory which were developed to study the classical Wigner matrices. Spectral properties 
of random adjacency matrix (II. 4p were examined in the limit n — ?■ oo both in numerical 
and theoretical physics studies. The first results on the spectral properties of sparse and 
dilute random matrices in the physical literature are related with the works [12) . |13| . |11] . 
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where equations for the limiting density of states of sparse random matrices were derived. In 
the papers and [5] a number of important results on the universality of the correlation 
functions and the Anderson localization transition were obtained. Unfortunately, these results 
were obtained with non rigorous replica and super symmetry methods. 

The first result on mathematical level of rigor for the matrices (jl.2p was obtained in [2] , 
where the eigenvalue distribution moments of the matrix (jl.2p with Pn = P were studied 
in the limit n — t- oo. It was shown that for any fixed natural m there exists nonrandom 
limiting moment lim„_>.oo '^~^Tr^'^ and these moments can be found from the system of 
certain recurrent relations. The results of |2j were generalized to the case of weighted random 
graphs in [8], where the resolvent of the adjacency matrix was studied and equations for the 
Stieltjes transform g{z) of the limiting eigenvalue distribution were derived rigorously (note, 
that the same equation for gaussian weights were obtained in [12], |13) . |11| by using the 
replica and the super symmetry approaches). But the limiting eigenvalue distribution, which 
is an analog of the low of large numbers of the probability theory, is only the first step in 
studies of linear eigenvalue statistics, corresponding to the test function ip 

Mn[ip]=Y,ip{Xi)=TTip{A). (1.3) 

Here and below {Aj}^^ are eigenvalues of the matrix A. The next step is to study the 
behavior of fluctuations of linear eigenvalue statistics. For the case of sparse random matrices 
this step was done in |17j with some modification of the method of 1^. It was shown in |17| 
that the random variable n^^/'^{J\fn[f] — E{Afn[f]}) converges in distribution to the gaussian 
random variable, as n — t- oo (here and below E{...} means the averaging with respect to all 

{o-ij}l<i<j<n) ■ 

The case of diluted matrices (p„ — )• oo) is less complicated technically than that with 
Pn = P- It was shown in [7J that in this case to have finite limits for E{A/'„[99]} one should 
consider the matrix A' = Aj ^/p^. Then it was proven in [7] that for integrable test functions 

2 

lim E{AA„[v3]} = ^ I (/p(A)v'4-A2dA, 

n,p„->-oo,p„/n-5>0 /TT J _2 

which coincides with the limits for the Wigner model [19]. Let us note that the method, 
used in [7], is rather similar to that for the Wigner model. But the problem to study the 
fluctuations of linear eigenvalue statistics usually is much more complicated than the problem 
to find the limiting eigenvalue distribution of random matrix ensemble. Even for the classical 
Wigner case the central limit theorem (CLT) for fluctuations of linear eigenvalue statistics 
was proven only recently in the series of papers with improving results |151 [H IH] 118] . 

In the present paper we prove CLT for fluctuations of linear eigenvalue statistics of diluted 
matrices, more precisely, we prove that the random variable (p„/n)^/2(A/'r,,[93] — E{A/'„[(/9]}) in 
the limit n^pn — ?• oo, pn/n — t- converges in distribution to the normal random variable. The 
method of the paper is a generalization of that of [18j. It allows us to prove CLT under rather 
weak assumptions on the test function ip (see Theorems [Tj and [2] below). 

It will be more convenient for us to study the matrix A = A/ y/p^ — Ej^/y/p^}, where 
E{...} means averaging with respect to all entries of A. It is easy to see that A differs from 
A! by the rank one matrix ^{A/ ^/p^}. So, everywhere below we will assume that the entries 
aij of A are distributed as 

^"^^ with probability Pn/n, 



aij = < 



/Vn n ^ 4^ 



n 



with probability 1 — pn/n, 
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Let us note that the case pn ~ an here corresponds to the Wigner ensemble, hence the model 
(jl.4p allows us to make "smooth transition" from the matrix studied in [8] to the Wigner 
matrix. 

Let us set our main notations. For any measurable function / we denote by E{/(^)} the 
averaging with respect to all random variables {ajj}i<j<j<„ and 

Var{/(^)} := E{|/(^) - E{f{A)}\'}. (L5) 

We denote also for any random variable ^ 

l = c = c-m}- 

Introduce the resolvent of A 

Gjk{z) = {A-z)-,^, Qzy^O, 7„(z) = TrG(z). (L6) 
In what follows it will be important for us that 

n 

\\G\\ < Yl l^^^f = (^^*)- ^ ll^ll' ^ l^^l"'' (1-7) 

i=i 

9(Ge,e)9z>0, VeGM". (L8) 

Here and everywhere below ||^|| means the operator norm of the matrix A. 

The main result of the paper is the central limit theorem for the linear eigenvalue statistics 
of any sufficiently smooth function ip which grows not faster than exponential at infinity. But 
we prove CLT first for the functions, which are smooth enough and decaying. Set 



Ml 



J{1 + 2\k\f'\$ik)\^dk, ${k) = e'''^vix)dx (1.9) 



and let Tig be the space of all function possessing the norm ||.||<j. 

Theorem 1 Consider the adjacency matrix l{1.4\ ) with pn — >• oo,p„/n — )• 0. Assume that the 
real valued function ip G Tig with s > 3/2 and that 

2 2-^2 



(^(/i)— =^d^ / 0. (1.10) 

o 

Then the random variable {pn/n)^''^J\fn[^] converges in distribution to a Gaussian random 
variable with zero mean and variance 

It is interesting to compare (|l.lip with that for the Wigner model 

M = n~'/''{wi,}l^^,, E{w^,} = 0, E{\w^j\^} = l, (i^j), 
we have (see [18]) 

1 /-^ /<^(Ai) -(^(A2)V (4-AiA2)dAidA2 



+ S ( / V^(M)^=^^/^ 1 + 



\ 2 / \ : 

2 2 - /i^ \ j f"^ (p{fi)fidfi\ 



27r2 \ ^'"^^ 71^7? ^1 47r2 \J^^ y^I^ I ' 
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Here K4 = n'^{E{M^j} - 3E^{Mfj}) = B{wfj} -3,W2 = nE{\Mii\'^}. One can see that (fTTTl) 
coincides with the term multiplying K4. This can be understood if we recall that in our case 

o 

= n'^(E{ajj} — 3E^{a?j}) ~ n/pn and we consider the random variable {pn/ny''^Afn[f], 

o 

while in the Wigner case one should consider J\fn [f] ■ 

One more interesting question is what is happening if the l.h.s. of (ll.lOp is zero. It is easy 

o 

to guess that in this case one have to change the normalization factor in front of J\fn [92] • But 
it could happen that the new expression for the limiting variance in this case will depend on 
the rate of convergence of Pn/n — >■ 0. We are going to study this situation in the future works. 

(c) 

Consider the set Tis of the functions, represented in the form 

93(A) = cosh(cA)^(A), ipeTis- (1.12) 



Theorem 2 Consider the adjacency matrix 1^1. 4\ ) with pn — ^ 00, Pn/n 0. Assume that the 
real valued function G Tii'^^ with some c > 0, s > 3/2 and II 1.10]) is satisfied. Then the 

o 

random variable {pn/n)^''^Mn[p>] converges in distribution to a Gaussian random variable with 
zero mean and variance U.ll\} . 



2 Proofs 

The proof follows the strategy developed in citeS:10 for the Wigner model. We start from the 
lemma 

Lemma 1 Let 'Jniz) be defined by U.6\) . Then for any z : > there exists a constant C 
such that 

^Var{7„(z)} < (^1^)' ^i^zt} < C/\^z\'\ (2.1) 

Moreover, for any e > we have 

^Var{7„(z)} < CE{|Gii|i+^}/|9z|3+% (2.2) 
n 

and for any smooth function F and any z : Qz > a 

n 

Var{n-i Vf(G,-,(z))| <n-i sup |F'(C)p. (2.3) 

j=l ^ C:0<3C,[C|<a-i 

Proof of Lemma[I\To prove (|2.ip we use the following proposition proven in f5] 

Proposition 1 Let ^a, a = l,...,i/ be independent random variables, assuming values in 
M™" and having probability laws Pa, a = 1, . . . ,1/ and let $ : R"^! x • • • x M™" C be a 
Borelian function. Set 

<J>a(ei,---,Ca) = j HCu---,Ca,Ca+l,...,Cu)Pa+l{d^a+l)...PAdCu) (2.4) 

SO that = ^, '^o = where E{...} denotes the expectation with respect to the 

product measure Pi ... P^. 

Then for any positive p > I there exists C , independent of u and such that 



V 



E{|$ - E{$}|2^'} < E{|$« - ^^-il^^'}. (2.5) 



a=l 
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Let $ = 7n(-z), = {0'aj}j<a- Denote also Ea{-} the averaging with respect to the random 
variables {aaj}j=i- Then it is easy to see that 



and by the Holder inequality 



Ea+l ■ ■ ■ 



Ell^^-^^^il^^'} = E{|i?„+i . . . Eni^-Eaim""} < = E{|$-Ei{$}|^f} 

(2.6) 

Define A^^^ (n — 1) X (n — 1) matrix which can be obtained from A if we remove from 
A the first line and the first column. Set also 



i=2 



We use the representations: 

Gn(z) = -(z + (G(i)a«,aW))-i 



(2.7) 



(2.8) 



(GWaW),(G«a«), 



z+(G«a(i),a(i)) 
Since G^^^ does not depend on a^^\ we have 

l + (G(i)G«a«,aW: 



. l + i?(z) , ^ (1 + Biz) 



z + (GWaW,a(i)) 



^(z) 



+^n^4(^/- 



(2.9) 



Hence, it suffices to estimate B{\B/A - Ei{S/A}|2} and E{\A-'^ - Ei{^-i}|2}. We show 
how to estimate the first expression. The second one can be estimated similarly. Denote by 
= C ~ for any random variable ^. Note that since for any a 



Ei{|C - a\'} = B^{\e,\'} + \a - Ei{e}|^ E.HCil'} < E^{\^ - a\'} 



(2.10) 



it suffices to estimate E{\B/A - Ei{S}/Ei{^}p} instead E{\B/A - Ei{B/A}\'^}. Then it 
is easy to see that 



B Ei{B} 



A Ei{A} 



Bl 



Al B 



< 



Bl 



Ei{A} 



+ 



Al 



•^zEiiA} 



Ei{^} Ei{^} A 
Here we used the relations that follow from the spectral theorem 

S>(GWa«,a«) =az(GWaW,G«aW), 9Tr G^^) = QzTV (G^G^^ 



(2.11) 



< \^z\-\ 



n 



-^Tr (GWgW* 



< \Qz\-\ 



|z + (G(i)a(i),a(i))| " ' |z + n-iTrG(i), 
The first relation yields, in particular, that |-B/^| < It is evident that 

Al = J2G\faua,j + J2G^\<^lir, 

i^j i 

Ei{|^?P} < G(p„n)-iTr(G(^)G(^)*). 



(2.12) 



(2.13) 
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In view of (I2.12P and (ll.7p we have 



^ I.. , „_irT...^(-ni. ^ ^PETTTTiJ- (2-14) 



|z + n-iTrG(i)| " |z + n-iTr " |Ei{yl}|^ 

Here in the first inequahty the numerator and the denominator D are just written as 
AT = iV^iVl-^i:» = D^D^-", then for {N / Df-'' the second inequality of (|2T^ is used, and 
then for the inequahty (ll.7p is used. Thus, in view of the second hne of (I2.13P 



E 



Ei{^} 



|z + n-iTrGW|2 " p„|Ei{^}| 



El 



Similarly 

Bl 



Ei{5} 



Tr(G(i)GWGW*GW*) p-^n-iTr (G^GW*) 

^ L/ ; — O " TTTTTr Cy - 



p„n|z + n-iTrG(i)|2 " |9z|2|z + n-iTY G(i) p " |Ei{A}|i+^ ' 

because, using the averaging with respect to {oij}, we obtain for the same bound as 

in the second line of (|2.13p . but with G*^^) replaced by (G^^))^. This gives the first inequality 
above. Then we use that Tr (G^^^G^G^i^G^^)*) < |9z|-2TV (G^^^G^^)*) (since HG^^^lP < 
|9z|-2 and finally use ([214]) . 

Then, the Jensen inequality |Ei{^}|~"^ < Ei{|A|^"'^}, and the relation = —Gii{z) 
yield 

Then (j23|) for p = 1 implies ([22]). Putting here e = we get (f2T]) . 

To prove the second inequality of (12. ip . we use (|2.5p for p = 2. In view of (j2.9p it is enough 
to check that 

E,{\Al\^} < Gp-2|9^|-4^ Ei{|i?°|n < Cp-^\Qz\-\ (2.15) 

The first relation here evidently follow from (j2.13p . if we take the fourth power and average 
with respect to {an}. The second one can be obtained similarly. 

To prove ([23]) we note first that ([23]) and ([22]) for $ = n"^ ^ FiGjj) yield 

Var{$} < nE{|$ - Ei{$}|^} < n'^E^^Yl i^i^jj) " ^i^^ff)) Q 

i 

<n-i sup |F'(C)Pe{(J;|G,,-G«|)\ 

if we take into account that ^ FiG^^j) does not depend on oij and hence may play the 
role of a in the inequality (|2.10p . Moreover, using (j2.8p and (j2.12p . we get 

^1 ^■■'■1- |z + GWaa),a(i))| ' ■ 

The above two bounds prove (|2.3p . □ 

Lemma [1] gives the bound for the variance of the linear eigenvalue statistics for the func- 
tions 99(A) = (A — z)~^. Now we are going to extend the bound to a wider class of test 
functions. For this aim we use Proposition below. We formulate it for the variance of linear 
eigenvalue statistics, but one can see easily that it can be applied also to a more general case 
even without reference to random matrix, see e.g. |14) . Proposition was proven in |18) . but 
for the completeness we give its proof here. We also would like to thank Prof. A.Soshnikov 
for the fruitful discussion on the proposition, which allows us to make proof the proof more 
simple. 
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Proposition 2 Let A be any random n x n matrix, A/'„[c/j] be its linear eigenvalue statistic 
U.3\) . and 7n(-z) be defined by U.6\) . Then 

Var{Mn[^]} < CsWipfJ dye-yy^'-^ Var{7„(a; + (2.16) 

JO J-oo 

where \\'^\\s is defined in \1.9^) . 

Remark 1 If the integral in the r.h.s. is equal to infinity, then the inequality is not interesting, 
hence we will assume that this integral is finite. 

Proof. Consider the operators Vg, V defined in tiie space of tlie Fourier transforms of tlie 
functions of the standard 



v.j{k) = {i+2\k\rf{k), 

Vf{k) = j dk'V{k,k')f{k'), V(fci,fc2) = Cov{W''i-^, W'^^^}. (2.17) 

It is easy to see that if we introduce the operator K := D^^VP^^ then 

Var{A/'„M} = (2^)-2(V^,(^) = (2^)-2(KP,(^, P,(^) 

< (2^)-2||i^|| . ||p,^||2 < (2^)-2||^||2TrK (2.18) 

Let us check that the operator K is indeed of the trace class in L2{^). Note first that in 
the Fourier space his kernel has the form 

k{ki,k2) = (l + 2|fcl|)-^V(A;l,fc2)(l + 2|fcl|)-^ 

with V{ki,k2) of (I2.17p . It is evident that X > in the operator sense, and K{ki,k2) is a 
continuous function of ki, /c2, since V(fci, /C2) can be written as a finite sum of the products of 
Fourier transforms of the positive unit measures, which are the distributions of eigenvalues of 
A. Moreover, we will prove below that 

K{k,k)dk< 00. (2.19) 

Then it follows from the inequality \K(ki,k2)^ < K{ki,ki)K{k2,k2) (which is valid for 
any continuous positive definite kernels) that K belongs to the Hilbert-Schmidt class, and 
therefore K has a basis {(pjik)}^^^, which is made from the continuous eigenfunctions with 
corresponding eigenvalues Xj > 0. Then if we consider a finite rank operator with the ker- 
nel K]\f{ki,k2) = 'Yl,^=i^j4'j{ki)<Pj{k2), we have K — Kj^ > in the operator sense, and 
K{ki,k2) — Kj^{ki,k2) is continuous, thus K{k,k) > K]\f{k,k) and 

N 



J2^j= [ KN{k,k)dk < I K{k, 
i=i 



k)dk. 



Since N is arbitrary and Xj > we obtain that K is a. trace class operator. 
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We are left to prove the inequality of (I2.19p . We have 
jK{k,k)dk = j{l + 2\k\)-^'V{k,k)dk 

= dye-yy^'-' J e-2|'=l?'V(fc, k)dk 

= dye- y / / dkidk2e'^''^~'''''>''V{ki,k2)e~^'''^y-\''^\y 

= ^ dye-yy^^-' J dxVar{AA„[P,(x - .)]} 
= ^ dye-yy^'-^ J dxVar{'^-fn{x + iy)}, 



where Py is the Poisson kernel 



Py(^) = ^tJtTITv (2-20) 



7r(x^ + y 
and we used that 

This relation combined with (j2.18p proves (j2.16p .D 

Now we are ready to prove the bound for the variance of linear eigenvalue statistics for a 
rather wide class of the test functions 



Lemma 2 If II95II3/2+Q: ^ cO; with any a > 0, th 



en 



^-^Y^viMnM] <CM\l,2+a (2.21) 

Proof. In view of Proposition [2] we need to estimate 

/oo 
Wa.T{-in{x + iy)}dx 
-00 

Take in (|2.2p e = a/2. Then we need to estimate 

/oo 
^{\Gii{x + iy)\'+^/^]dx. 
-00 

Use the spectral representation 

Gu=[^^^^, were Arn(A) = f; |4'^I^1a(A.) 
J X-x-iy f-' 

with = (V'f^.-.jVn^) being an eigenvector of A, corresponding the eigenvalue A^, i.e. 
Aip^^^ = Xkip^^\ Then the Jensen inequality with respect to Nii{dX) yields 

II i^i" ^ -^"^ ^ lyc (I. - /:';;aw./. <- ^i"!-"' 

Taking s = 3/2 + a in (j2.16p we get 

^Var{AA„M} < H e-yy'+^'^y-'-^dy < C\Ml,_^^. 

□ 

The next lemma is technical one. We accumulate relations which we need to prove CLT. 
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Lemma 3 Using notations of Ii2.9\) we have uniformly in zi,Z2 ■ 9-Zi,2 > ^ with any a > 0; 

B,{\Ar} = 0{p-% Ei{|S°|n = 0(p-2), (2.22) 

(Ei{A-n)° = - (l + 0{p-') + Oipjn)) -^^^ + r, (2.23) 
with S{|r°|2} < C/n^ + C/pln, 

p.Ei{^°(zi)^°(z2)} = ^E4'^(^i)gSP(^2) +Pnt\zi)t\z2)/n', (2.24) 

i 

PnBi{A°izi)B°iz2)}=Pn:^E,{A°{z,)A°{z2)}, (2.25) 

az2 

VarKEi{A°(zi)A°(z2)}} = 0{n-^), Var{p„Ei{A°(zi)i?°(z2)}} = O(n-i), (2.26) 

B{\7n\z)-Uz)\'} = 0{p-'). (2.27) 



Var{Gy)(zi)} = 0{p-'), \B{g\1\z,)} - E{G,,(zi)}| = 0{p-'), (2.28) 

|E{7i')(z)}/n - f{z)\ = 0{n-'), \E-'{A{z)} + f{z)\ = 0{p-'), (2.29) 

where 

f{z) = \{^/^A-z). (2.30) 

Proof. Note that since ^^^(G^^^m, m) > 0, we can use the bound 

13=^1 > ^ l^^^l < < (2.31) 
Relations (12.221) . (12.241) . and (12.251) follow from the representation 

^° = E + E 4^(4)° + n-Si'\z) = Al + n-'t\^), (2.32) 

and Lemma □ (see ([211) and ([2J5|) ). 

The first bound of (|2.26p follows from (I2.24p and (I2.3p for F{z) = . The second bound of 
()2.26p follows from (|2.25p and the first relation of (|2.26p , if we use the fact that the variance of 
the derivative of an analytic function by the Cauchy theorem can be bounded by the variance 
of the initial function. 

Relations (I2.27P follow from the representation (see (I2.8P ) 



7r(z)-7„(^) = (^-')° + (i?^-V 
and (|2.34p . The first relation of (j2.28p is the analog of the relation 

Var{G,,(zi)} = Var{Gii(zi)} = 0{j>-^) (2.33) 
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if in the latter we replace the matrix A by A^^^ . But since Gii(zi) = —A~^{zi), (I2.33P 
follows from (j2.22p and (j2.3ip . The second relation of (|2.28p follows from the symmetry of 
the problem and (I2.9P 



B{G^\z,)} - B{Gu{zi)} = ^nin - 1^^ - Gn} = -l^E{i?M} = 0{n-'). 

The first relation of (j2.29p follows from the above bound for n~^E{7„ — 7^^^} and the estimate 
(see IZl) 

n-iE{7„}-/(z) = 0(p-i). 
The second relation of (I2.29P is the corollary of the above estimate and the representation 

A-i = E-i{A} - ^°E-2{^} + (A°)2A-^E-2{^}, (2.34) 

which implies 

E{A(z)}-i = E{^(z)-i} + 0{\^{A{z)]) = -E{Gn} + 0(p-^) = -n-iE{7„} + 0{p-'). 
We are left to prove ^a?M . Set 

A = z + Y,G^^al. (2.35) 
Using the analog of (I2.34P for A and A, we write first 

A-^ = A-^ - A-^{A - I) + ri, n = A'^A~\A - Af. 

We have 

A-A = Y^ G\fa,,a,„ B{\n\^} < jJ^B{\A - A\^} < Cjr? + Cjuvl. 

■_L- 

Moreover, the analog of (|2.34p for A yields 

B^{A-\A - A)} = Bl\A}Bi{{A - A)} 

- 2Er'{A}Ei{(^ - I)(I - Ei{A})} + rs = r2 
r2 = Ei{Er2{A}(l-2 ^ 2B^\A]A-^) [a - l) (l- Ei{^})'}. 

Since A - Ei{A} = ^ G^Vh " Ei{a?J) we have 

E{ri} < CE{|^ - - Ei{A}|4} < C/np^ . 

Hence we have proved that 

Ei{A-i} = Ei{I-i}+f, < C/n2 + C7/np2. (2.36) 

Using (II. 4p we can write 

poo /-^\ 

Jo 

= dte'^V'P-'^-' 1^" E (e**^-'(i/^"-2/") - 1) } + O(n-i) 

= j dte''V^- (i-P"/")/"exp|m-i^F(G£\t)| + 0(n-i), 
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where 

F{x,t) = p„(^e^*-(i/Pn-2/n) _ 1 _ it^^i/p^ _ 2/n)) , 
Then in view of (|2.3p . since supQ^^g \F^{x,t)\ < C\t\p~^, we obtain 

.E.{I-} ^ f ..e-e-^'(-^-^/")/"exp{..-5:E{n4),,}} 
.(l + 0[n-'^F%Gil\t)))+0{n'') 

= J (ite^^*e^*^" '(i-P"/")/" exp jm"! ^ E{F(g[^\ t)}} + r', 
E{|rf} <C/n2 + C/np2. 

Finally, replacing similarly to the above 7^^'' by E{7^^''} in the exponent, we get 

iBi{A~^} = y"°°dte^^*e**^^T"'>(i-^'"/")/"exp{m"i^E{F(G2\t)}} 
(1 + itn-\ji^^ni - pjn) + 0((n-i(7(i))°)2) + r'. 

Taking (Ei{^^^})°, we can see that the term which corresponds to 1 in the r.h.s. disappears, 
and since 'E{F{Gii,t)} = 0{p-^), the coefficient in front of (7!^^ equals 

' i^^^zt^im{,i!^}(l-p„/n)/n ^ E{F{G^^ , t)}}dt 





= i{z + E{^W}/nr\l + 0{p-') + 0{pjn)) 
= -zE-2{A}(l + 0{p-') + Oipjn)). 

In view of dHJ and (I236D we obtain 
□ 

Proof of Theorem [7J We prove first Theorem [1] for the function ip^^ of the form 



Pr,*ipo, j \M>^)\dX<C <oo, (2.37) 

where iP^ is the Poisson kernel (see (j2.20p ) and ipQ is a real valued function from Li(M). One 
can see easily that 

M°[vr,] = {Tnpr^iA))" = <^o(/u)9(TrG(/i + ir?)) V 

= ^ / ¥'o(/^)(7n(^M) - ln{zti))dl^, z^ = ix + if]. (2.38) 



Set 



dn=f-V^', Z„(x) = Eje^^^-I^-l/'in}^ (2.39) 



e(x) = e^^-^-I'^"]/'^-, y„(z,x) = d-^E{TrG(z)e°(x)}. 



Then it is easy to see that 



-^Znix) = I foifJ-)iYn{Zf,,x) - Yn{z ^, x))d^i. (2.40) 
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On the other hand, using the symmetry of the problem and the notations of (12. 9p . we have 

Yn{z,x) = d-^B{TvG{z)e°{x)} = nd;;^B{Gu{z)e°{x)} (2.41) 
= -nd~'E{{A-^yei{x)} - nd-'E{{A-^y {e{x) - ei{x))} = Ti + T^, 

where 

e,{x) = e-(A/'^-\[^.])V'^., (AA«iK])° = (Trv^.^^^)))" = ^ | df^Mf^)^ln (z,)- 

Since ei(x) does not depend on {an}, using that E{...} = E{Ei{...}}, we obtain in view of 
the above representation and (j2.23p 

Ti = d-^EMHz)yei{x)}/E'{A}[l + 0{d-') + 0{p-')) + 0{d-'). 

Write 

eix) - ei(x) = ^ 1 ^o(/^)(9(7n " 7n^'^)ei(x) + ©((t^ - 7n^'V)^^/dn)^i/.. (2.42) 

Then ^23), the relations |e(x)| = |ei(rE)| = 1, and yield 

d~'\BM^ye,{x)} - E{^:e{x)}\ < Cd~'E{\{j^^^y - 7^1(1 + \x\K\d~')} 

Hence we obtain 

T, = Yn{z,x)/E'{A}(l + 0{d-') + 0{p-')) +0{d-^p^^/^). (2.43) 

To compute T2 we use (I2.42p . Then, taking into account (j2.27p . we conclude that the term 
0(rid~'^(7° — (7!^'*)°)^) gives the contribution 0{d~^). Then, since 61(2;) does not depend on 
{aii}, we average first with respect to {an} and obtain in view of (12.91) 

T2 = -'^l df,M^^)ES^{A-'yiz)ei{x)-s{j:{z,) - {7^r^Hz,)y)} + 0{d-') 

= I d^^o(/x)E|ei(x)Ei{(A-i)°(z)9(7:(z^)-(7i^)(^M))°)}}+0(d-i) 

= <i^^o(/i)E|ei(x)Ei{(^-i)°(z)9((l + S(z^))A-i(z^))°}|+0(d-i). 

Using (I2.34P and (I2.22p . we conclude that only linear terms with respect to B° and A° give 
non vanishing contribution, hence we obtain 

Dniz,z^) := pnE^{{A-'y{z)(^il + B{z^))A-\z^)y] 

= PnE-^{A{z)}E-^A{z^)}(^l + E{B{z^)})Ei{A^z)A^z^)} 

-PnE-^Aiz)}E~'{A{z,)}E^{A%z)B%z,)} + Oip~'/') 
= f{z)fiz,){l + f'{z,)) + f\z)f{z,)f{z,) + 0{p-y'). 

Here we used first ^(TI^ and (12:251) to express Ei{A°{z)A°{z^)} and Ei{A°{z)B°{z^)} in 
terms oiG^^{z), and -^G^^{z^), and then[228]) combined with (l2:29]l to replace Ei{A°{z)A°{z, 
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by f{z)f{Zf,) and Bi{A° {z)B° (z^)} by f{z)f{z^). Moreover, we used dOO]) to replace 
E-i{y4(z)} by -f{z). Hence 

D^{z, z,) = [f{z)f{z,){l + f'{z,)) + f{z)f(z,)f'{z,)) + 0{p-'/'). (2.44) 

In addition, similarly to (I2.42P we have 

E{ei(x)} = Z„(x) + 0(d-i). 
Hence, relations (|2.4ip - (|2.44p imply 

Yn{z,x) = f'^{z)Yn{z,x) + ixZn{x) I dflifioifJ')- 



p2/^^^z ™\ I l\i ^ ,,\-^n{z, z^) Dn{z,z^) 



2i-K 



+ 0(p-V2) + o(d-i^ 



y.(.,x) = ^xZr.ix) j d;.v^oM ^"^"'""\./"^^'""^ + 0(P.^/') + OK'), (2.45) 



Cn{z, Zfj) : 
Using that 



l-/2(z)' 



/(^)(m + i) = :r^ = -^, /'- ^^^^ 



we can transform Cn{z, z^) to the form 

C„(2,2„) = Clz.z^) + Ofc"2) + 0{d-"^) (246) 



(^2 _ 4)1/2(^2 _ 4)1/2- 

Taking into account (j2.40p . (|2.45p . and (|2.46p . we obtain the equation 

^Zn{x) = -xV[^^]Zn{x) + 0{p-^'^) + 0((i-i) (2.47) 



+ C{'zj;;,zj^fjdfiidii2- 

Formulas (121^61) and (j2:i71) imply that 



^[^'^] = ^(/ ^°(^)^ (/i 4)1/2 ^^) =i(/ ^°(^)^( (/-4)V2 -"^)'^^) 
where we used also the well known relations 



1 f'^ dX 1 1 /"^ dA 



2 (z - A)V4^ (^2 _ 4)1/2 ' ^ ^/4^ 



1. 



Now if we consider 

Z„(x) = e^'^['^''l/2z„(x), 
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then (I2.47P yields that for any |x| < C 

and since Z„(0) = Zn{0) = 1, we obtain uniformly in x < C 

Zn{x) = l + Oip7j/') + 0{d-') 
^Zn{x) = e-'^K]/2 + 0(p-i/2) + 0(^-1), (2.48) 

Thus, we have proved CLT for the functions of the form (j2.37p . To extend CLT to a wider 
class of functions we use 

n 

Proposition 3 Let {£,[^^^^^1 be a triangular array of random variables, Mn['^] = ^('^j'"^) 

1=1 

be its linear statistics, corresponding to a test function 99 : M — ?■ M, and 

VnM = Var{d-^MnM} 

be the variance of J\fn[(p], where {dn}^=i is some bounded from above sequence of numbers. 
Assume that 

(a) there exists a vector space C endowed with a norm ||...|| and such that Vn is defined 
on C and admits the bound 

Vn[v]<C\M^ yipeC, (2.49) 

where C does not depend on n; 

(b) there exists a dense linear manifold £1 C £ such that the Central Limit Theorem is 

valid for Mn[^], ^ G £1, i-e., if Zn[xip] = Eje"''^" A^nMj ig ^/^g characteristic function of 
— 1/2 ° 

dn Afnl'p], then there exists a continuous quadratic functional V : £1 — )• M+ such that we 
have uniformly in x, varying on any compact interval 

lim Zn[xip] = e-^'^[^l/2, V99 G £1; (2.50) 
Then V admits a continuous extension to £ and Central Limit Theorem is valid for allj\fn[<f], 

Proof. Let {fk} be a sequence of elements of £1 converging to ip £ C. We have then in 

o 

view of the inequality |e*" — e**| < |a — 6|, the linearity of Mn[^] in ^, the Schwarz inequality, 
and (loo]) : 



Zn{x(p) - Zn{x(p)\ 



<fi='Pk 



< \x\E 



djMn['p]-d,,^Mn['Pk] 



(2.51) 



< \x\Yar'^/^{djAfn[ip - y^k]} < C\x\ \\^ - (pk\\- 



Now, passing first to the limit n — )• 00 and then /c — )• 00, we obtain the assertion of the 
proposition. □ 

Let us show now that hypothesis (a) and (b) of Proposition [3] are fulfilled in some vector 
space. Consider the space Tig of all functions with the norm (II. 9p and set £ = Tig D Li (M) 
and 

'|<^(A)|dA + ||</.|U = ||<^||i + (2.52) 
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Then for s > 3/2 Lemma [2] guarantees that assumption (a) of Proposition |3] is fulfilled. 
Moreover, the Lebesgue theorem about the dominated convergence yields that 

\\^-^*P^\\l<C j \l-e-'^\^\\^{l + 2\k\f'\^{k)\'^dk^Q, 7]^0. 

Hence the set of the functions (p * Pr^ is dense in C with respect to the norm ||.||s. Thus, if 
we prove that the set of the functions (p* Pr/ is dense in C with respect to the norm 1 1 • 1 1 (k) j 
then (I2.48P will imply assumption (b) of Proposition [3l 

It is easy to see that the set of all functions with finite supports, possessing the norm 
(j2.52p . is dense in C with respect to this norm. Hence we need only to prove that if (/? S Tig 
and has a finite support [—A, A], then 



J \ip{X) - if* Pr^{X)\dX ^0, rj^O. 



But 

[ \ip{X)-ip*P^{X)\dX= ( [ +[ ]\^{X)-^*Pr,{X)\dX = h+l2 

J \J|A|<A+1 J\\\>A+1/ 



We have for I2 



J\\\>A+i J\x\<a{>^- fj-r + v 



and for Ii we use the inequalities: 



Ii<2|^ + l| sup \(p{X)-^*Pr,{X)\<C\A + l\ [ \l-e-'^\''\\\lp{k)\dk. (2.54) 
lAKA+l J 



But since 



(I2.54P the Lebesgue theorem on the dominated convergence implies that /i — )• 0, as 77 — >■ 0. 
Combining this with (|2.54p we get that the set of all functions with finite supports, possessing 
the norm (I2.52P is dense in £. As it was mentioned above this implies that the set of the 
functions ip * P^ with € £ is dense in C with respect to the norm ||-||l^(ir) and in view of 
(j2.48p proves assumption (b) of the proposition. □ 

Proof of Theorem [H Let us note first that in the case when 

Pn > C^:log^^^ n, n — ;> 00, (2.55) 

the proof of Theorem [2] is rather simple. By the method of [3j one can prove the estimate 

3 

n-^EiTiA^""} < C^il + — )n-^E{Tr^2m-2|_ 

Pn 

Then under condition (|2.55p it is easy to get the bound, valid for sufficiently big K: 

Prob{||^|| >K}< inf nE{Tr (^/i^)2'"} < e^p{-pl/^ log{K/2C) + logn} ^ 0,n ^ 00. 



15 



Then, for any (/? G 'Hi'^\ if we consider a smooth function ip^^'^ £ Tig with a finite support and 
such that ip^^^X) = (^(A), |A| < K, then evidently 

|j.|gixM?M/rfn} _E{e^^^n[<^*'"]M^}| < Prob{||^||K} 

+ d-i|E{AA„[vp]} - E{AA„[(^(^')]}| ^0, n ^ oo. 

Thus, we can derive Theorem [2] from Theorem [T] almost immediately. 

But if the inequality (12.55^ is not fulfilled, then the proof of Theorem[2]is more complicated. 
It is based on the bound which is the analog of (|2.16p 



^Yar{KM} < C{c)M\s, 'P e ni'\ (2.56) 
n 



where (p{X) = </'(A) cosh ^{cX). The main step here is the lemma, which is the generalization 
of Lemma [T] 

Lemma 4 Denote by jn^ = Tr G(2;)e'^'^. Then for any 1 > e > 

^Var{#} < C(c,e)E{|Gii|i+^}/|9z|3+^ (2.57) 
n 

Proof. According to Proposition [1] it is enough to prove that 

E{|7(^) -Ei{7(^)}|2} < C{c,e)E{\G^,\'+'}/\-szf+'pn. (2.58) 

Let us set 

G(i)(z) = (^(1) - zr\ 7(1^) = TrG(i)(z)e^-4'''. 

Note that differently from the proof of Proposition [1] here and below we denote by A^^^ the 
n X n matrix whose first line and column are zero and the other entries coincide with those 
of A. We also denote a^^^ = (0, ai2, . . . , ain)- Then we can write 

7(^) - Ei{7('=)} =7^^) - 7^'^) - Ei{7(^) - it^}, 

7W -7(1-) =Tr(G(z) -G«(z))e^^*'^ +TrG«(z)(e^-^-e^^*'') 
+ Tr (G(z) - GW(z))(e"-^ - e"-^*'^) = / + // + m. 

Let us use the formulas 

(G(z) - G«(z))n = - A-\ {G(z) - G«(z))i. = -^-i(G«a«)„ 
(G(z) - G«(z)),, = -A-i(G«a«),(G«a«),-, i,j > 2, 



gC^_gC^W)^_^2 I" *d^(e^*^«a«),(e^^-^'''G«),(e^(i-*-^)-^)n, i,j>2. 



lo Jo 



Here A is defined in (j2.9p and to obtain the last two lines we have used the Duhamel formula, 
valid for any matrices M and M^^^: 



Jo 

Moreover, we have taken into account that 

A^=A^ = G^ = ie^^^'^)a=0, 4?=0, gS = -.-i, (e*-^'^^)n = L 
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Hence we have 

II = - z-\{e'^U - 1) + c2 f dt dr{e<^~'-^^^)u{G^^^e''^'"\^^\e'^^''\^'-^) 

Jo Jo 

Jo 



111 



,2 A~l 



c'A 



1 i-i-t 

dt 
Jo 
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Thus, denoting and using the Schwarz inequahty and (j2.10p - 

(f2TT]l . we get for / 



E{|/ - Ei{/}|2} < 3B{\A-' - Ei{A-i}|2} + 3e| ^^^Q 



+ 3E 



B^'^'^Al 



ABi{A} 



Averaging with respect to {an} and then using the Holder inequahty, we get 



El 



B 



(c)o 



Ei{A} 



< C- 



-iTr|G«|4e 



Pn|Ei{A}P 



< c- 



-iTr|G«|2e 



p„|9z|2|Ei{A}|^ 



— TTTTT^ Z- ■ ^ . ^ 



p„|9z|2|Ei{A}P 



p„|9z|3+^|Ei{A}|i+= 



where \G^^'^f = G^^^*G^^\ Similarly, using that QA = + {\G^^^\'^a^^\a^^^)) (see (f2J2D l. 
we obtain 



El 



^Ei{A} 



< El 



f / |2a(i), Q(i))i-^(e2'^-4^^V^|G(i) |2fl(i), a(i))^| A° 



|(l + (|GW|2a(i),a(i)))|Ei{^}| 



^ Ei{(e2^-^'^V^ad),a(i))2^|Af|2} ^ E}-^^{|A^ |2}Ef {(e^^-^'^V^gW, aW)|A^|2} 



|9z|2|Ei{^}|2 



|9z|2|Ei{A}p 



(n-iTV |G(i) |2)i-2-(n-i| |G(i) | pTr e2-^<'V^)2e (^-ixr e2c^(i)/.)2. 



Pn|9^P|Ei{^}|2 ^ p„|9z|3+2e|Ei{^}|l+2£' 

The terms with // and /// can be estimated similarly, if we use also the bound 

To prove (j2.59p . we prove first that 

Ei{((^'")ii - Ei{(^™)n})'} < Gm22-E}/'{(^4— 2)n}/p„. 
It is easy to see that 

[m/2] 

k=l iiH Hk='m.-2k 



(2.59) 



(2.60) 



(2.61) 
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with and a^^^ of (HZ]). Thus, using that for aU / E{a\j} > 0, we have 
Eij (s('i) . . . SC*) - Ei{S('i)} . . . Ei{S(^fc)})^} 

<A;^Ej/2{(s(^^)-Ei{S(^^)})'}E}/2{n(S('»))'} 



k 

< 



Taking the sum as in (j2.6ip and using the Schwarz inequahty, we obtain (j2.60p . The Taylor 
expansion, the Schwarz inequahty, and (I2.60p imply (12.59P : 



{{e^^)u - Ei{(e^^)n})n < Pn' E T^^^l (('^'")" " Eill-^'")!!}) } 

m=l 

m=l ^ m=l ' 



4m-2\ 



Lemma m is proven. □ 

The next step is the analog of Proposition [2] 

Proposition 4 For any 99 G Us 



Var{AA„M} <C,||(p(^)||2 / dye-^y'''-^ Var{ji^\x + iy)}dx, (2.62) 

JO J-00 
where <^(A) = (/3(A) cosh~"'^(cA) and is defined in i[1.9\) . 

The proof of Proposition |4] coincides with that of Proposition [21 if we replace the operator V 
of (I2.18P by the operator V'^'^-' whose Fourier transform has the kernel 

W){ki,k2) = Cov{TYcosh(c^)e*'=i-^,Trcosh(c^)e*'=2^}, Var{AA„[(^]} = (V^'^)^,^). 

Now one can derive (j2.56p from Lemma [J] and Proposition 2] by the same argument that 
we used in Lemma [2] to derive (I2.2ip from Lemma [T] and Proposition [2j 

Having in mind the bound (I2.56P , we can derive Theorem [2] from Proposition [3l if we are 

(c) 

able to prove CLT for some dense subset of Tig , e.g., for ip with finite supports. But if 99 has 

(c) 

a finite support and belongs to Tig , it belongs also automatically to Tis, thus we can apply 
Theorem [1] to it. This completes the proof of Theorem [2l 
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